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Abstract. We present classical molecular dynamics calculations of the behavior of 
copper and gold nanoparticles on a graphene sheet, sheared with constant applied force 
Fa- The force Fg acting on the particle from the substrate depends on the material 
of the nanoparticles (Au or Cu), and exhibits a sawtooth dependency on time, which 
we attribute to local commensurability between the metal nanoparticle surface atomic 
positions with the graphene lattice. The time-averaged value of (the friction force) 
acting on Au nanoparticles increases linearly with contact area, with slopes close to 
the experimentally observable ones. A qualitative model is proposed to explain the 
observed results. 
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1. Introduction 

The atomic force microscope (AFM) and surface force apparatus have greatly expanded 
the possibihties to investigate tribological properties of nanoscale objects In 
particular, the AFM has been used to study the frictional interaction between 
nanoparticles (NPs) on atomically flat substrates. These represent very well deflned 
experiments where both the size and shape of the nanoparticles and the nature of 
the substrate can be studied with atomic resolution both before and after sliding the 
particles |2|. 

A series of experiments concerned with the shear of NPs have been carried out 
recently [TJ[2]j3]|4][5 . The measurements suggest that the friction force acting on antimony 



NPs deposited on highly oriented pyrolytic graphite (HOPG) depends linearly on the 
contact area of the NP. Another interesting result is the discovery that some nanoislands 
with the contact area up to 9 • 10^ nm^ show very low friction forces, while other similar 
NPs exhibit much larger friction forces. 

In spite of some theoretical eflForts presented in literature [2','6], no deflnitive 
explanation of the experimental results has been proposed so far. Classical molecular 
dynamics (MD) simulations of silver and nickel NPs adsorbed on graphene and sheared 
with constant external force Fa were presented in Ref. f?^, where it was shown that 
the friction force acting on the nanoislands changes linearly with contact area. The 
friction force is quite sensitive to the material the NP is made of. In particular, the 
results suggest that Ni NPs on graphene experience higher friction than the Ag NPs. 
Moreover, a sawtooth form of the substrate shear force Fg acting on the nanoparticle, 
versus the lateral component of the center of mass position, was found for Ni while the 
corresponding curves for Ag NPs turned out to be rather irregular. Such a behavior 
was attributed to the local commensurability of the counterface of Ni NPs due to the 
close match of the nearest-neighbor distance in metal and the lattice constant of the 
graphene substrate. 

In Ref. [TJ only two metals were considered, and additional investigations are 
required to test how general the previous computational results are, and to prepare 
the foundations for a quantitative theory. With this aim we perform classical MD 
simulations described in this paper. We pursue two main goals. First, we want to assure 
the results related to the influence of the material of NPs on their tribological behavior. 
It is clear that short-range order of atoms can play a signiflcant role in determining 
the friction force |7|. If the distance a between nearest-neighbors in metal is close to 
the graphene lattice constant of 0.246 nm |8|, then local commensurability of atoms is 
possible, leading to a sawtooth type dependence of the force Fg acting on the NP from 
the substrate (graphene). The time-average (Fg(t)) = Ff of the (lateral) substrate force 
is the friction force. The commensurability can happen for Ni where a = 0.249 nm |9|. 
On the other hand, for metals with larger values of a commensurability should not occur 
and this was observed for Ag with a = 0.289 nm |7|. To conflrm this hypothesis we 
consider here two other metals: Cu with a = 0.256 nm and Au with a = 0.289 nm [91. 
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The paper is organized as follows. Section [2] describes the atomistic model. Results 
of the calculations and the discussion are given in section |3| Brief conclusions are 
presented in section [4| 



2. Model 

A detailed description of the model can be found in Ref. Here we provide only 
a general sketch and the details which are specific to the current calculations. The 
substrate is a graphene sheet parallel to the x^-plane, with zigzag and armchair edges 
along the x- and ^-directions, respectively (see figure [l| all snapshots in this work were 
produced with VMD software \10^). Boundary carbon atoms along the perimeter of 
the graphene layer are held rigid in order to provide the fixed position of the layer. 
Copper and gold nanoislands containing from 5000 up to 30000 atoms are considered. 
Depending on the size of the nanoparticle, the x x y dimensions of the graphene sheet 
varied from about 19.68 x 17.04 nm to 36.40 x 31.52 nm. The total number of atoms in 
the calculations varies from 17800 to 73808. 




Figure 1. Snapshots of the formed nanoparticle containing 22000 Cu atoms: 
perspective (a) and top (b) views. 



Interactions between carbon atoms in graphene are described by the harmonic 
potential 11 . Embedded atom method (EAM) potential is employed for forces between 
metal atoms 12 . For the metal-carbon interaction the pairwise 6-12 Lennard- Jones (LJ) 
potential is used 11 . The simulation code is implemented using the NVIDIA CUDA 



platform 13, 14] which allowed us to carry out the computations on a single graphics 
processing unit (GPU) NVIDIA GeForce GTX 460. Algorithms for GPU based on the 
neighborlist technique (l5 , 16 from Ref. (l4| with our own algorithm for binning atoms 
into cells are employed. The equations of motion are integrated using the leapfrog 
method fiE, 16 with a time step At = 0.2 fs. 



3. Results and Discussion 



NPs are obtained using the procedure similar to the dewetting of thin metallic films by 
thermal treatments 17 and described in Ref. |7|. Briefly, a rectangular slab of metal 



Atomistic modeling of friction of Cu and Au nanoparticles on a graphene sheet 



4 



atoms is placed on the graphene sheet and the simulation starts. The metal atoms 
assemble into a more compact structure, which is more energetically favorable than the 
rectangular slab. At an appropriate time moment the system is cooled down in order 
to get a hemispherical shape of the NP and to prevent it from forming a ball. After 
the equilibrium has been reached, the external force Fa is applied along the x-direction 
(zigzag edge of graphene), and the NP starts to move in this direction. This force is 
distributed uniformly on all metal atoms with the x < Xqm^ where Xqm is the NP 
center of mass position along the x-direction. We perform several series of calculations 
with different values of the applied force Fa acting on the NPs of the same size. Where 
it is not stated explicitly, the Fa acting on Cu NPs consisting of 5000 and 29000 atoms is 
equal to 2.31 nN and 15.29 nN, respectively. For Au NP consisting of 10000 and 25000 
atoms Fa is 5.07 nN and 13.05 nN, respectively. We note that in our MD calculations the 
NPs are accelerating since in general the applied force Fa is not equal to the substrate 
force Fg. 




time (0.2 fs) 



Figure 2. Time dependencies of the system temperature T, lateral position Xcm and 
velocity Vx of the center of mass of the NP, total applied force Fa, substrate force Fg, 
and lateral dimensions Lx and Ly obtained for Cu nanoisland containing 19000 atoms. 



Figure [2] shows the time dependency of the quantities measured during the 
simulation of a copper NP containing 19000 atoms. Note that once the motion of 
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the NP is initiated, the values of Vx and Xqm increase nearly linearly and quadrat ically 
with time, respectively. This is possible only if the friction force acting on the NP is 
independent of the sliding velocity (in the studied velocity range) so that and Fa — Ff is 
time independent. However, the substrate force Fs{t) acting on the NP is not constant, 
but exhibits a sawtooth shape with increasing frequency (in time) as the particle moves 
faster, which indicates an atomistic stick-slip motion of the NP. We note that using 
MD simulations (or any other exact numerical simulation method) one cannot study 
the low NP sliding velocities (of order ^ /xm/s) prevailing in the AFM measurements, 
but the numerical results presented below indicate a very weak velocity dependency in 
the studied velocity interval. 




(0.1 nm) 



(0.1 nm) 



Figure 3. Substrate force versus the lateral position of the center of mass of the 
nanoparticles: Cu with 5000 (a) and 29000 atoms (b), Au with 10000 (c) and 25000 
(d) 



To examine the sliding dynamics in more detail, in Fig. [3] we show Fg vs Xqm for 
several Cu and Au NPs. While for the majority of Cu NPs a sawtooth dependence of 
the substrate force is observed and Fg peaks are more or less regular, the Fg vs Xqm 
curves for Au have more irregular form. The distances between the peaks for Cu are 
close to the value of the lattice constant of graphene. Note that Fg takes both positive 
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and negative values. Thus the force acting on the NP from the graphene changes its 
direction during the shding of the NP. On the average, the amphtude of the peaks for Cu 
is larger than for Au. The substrate force for the Cu nanoparticle exhibits fluctuations 
with a periodicity close to the lattice constant as expected for a nearly commensurate 
contact. For the Au nanoparticles the substrate force fluctuates much more rapidly with 
the displacement coordinate Xcm, which is expected for a nearly incommensurate (or 
high-order commensurate) interface [18j . 




^cm(0.1 nm) Xcm(0.1 nm) 



Figure 4. Substrate force versus the lateral position of the center of mass of the Au 
nanoislands at low applied force: 10000 atoms (Fa = 0.74 nN) (a) and 25000 atoms 
(Fa = 2.08 nN) (b) 




Figure 5. Friction force versus contact area of Au nanoparticles for different applied 
forces and hence different sliding velocities. 

Figure [4] presents the dependence of the substrate force Fg on the center of mass 
coordinate, for two different Au NPs with different applied forces Fa. Note that the 
curves exhibit "flne structure" for all sizes of the NPs, i.e. the peaks are split into 
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high-frequency subpeaks, and that the stick-shp osciUations are reduced or disappear 
for the largest Xcm corresponding to the largest sliding velocities. To explain such a 
behavior one can assume that only at low enough speeds the atoms on the surface of NPs 
can adjust to the potential energy landscape of the graphene layer. In this way local 
atomic commensurability patches can form at the interface, resulting in atomic-scale 
oscillations in the substrate force acting on the NP. In contrast, at high sliding speed 
the interfacial atoms of the NP do not have enough time to adjust to the corrugated 
substrate potential, and no local commensurate regions or patches can form. This can 
be the reason for the irregular form of the substrate force dependencies for Au NPs at 
high sliding velocities. 

Figure [5] shows the dependence of the friction force Ff = (Fg) on the contact area of 
the NPs. It is noteworthy that the change in applied force does not significantly alter the 
average value of the friction force, i.e., the friction force is nearly velocity independent in 
the studied velocity interval. The slopes of the linear fits to the friction force data for Au 
NPs are 2.29 and 2.03 pN/nm^ at higher and lower load, respectively. They exceed the 
values of 1.21 pN/nm^ obtained in MD simulations for silver |7| and 1.04 pN/nm^ from 
the experiments (3j|^. However, the experiments are performed at much lower sliding 
velocities than in the MD simulations, and the fact that the slope for the lower applied 
force is slightly smaller than for the higher force is consistent with a smaller slope at 
the AFM sliding velocities. Poor statistics did not allow us to reliably determine the 
slopes for copper NPs. We note that if the NP atomic positions at the interface would 
be strictly incommensurate with the graphene lattice one would expect a friction force 
roughly proportional to the linear size of the NP contact area. The fact that the friction 
force is proportional to the contact area instead is consistent with a model where in the 
contact region the metal atoms form commensurate patches separated by domain walls. 

The structure of NPs can be analyzed using radial distribution function (RDF) 
measured at different time moments, and from snapshots of the contact surface. Typical 
RDFs are shown in figure [6| As can be seen, RDF has a completely blurred shape after 
the formation of the NP for both metals, indicating a disordered structure of the NPs. 
After the cooling phase some ordering occurs. It is manifested in the formation of 
the first peak, corresponding to the distance between the nearest neighbors of about 
0.256 nm and 0.289 nm for Cu and Au, respectively. Additional smaller peaks also 
appear in the figure. However, they are very diffusive compared to the ideal bulk crystal. 
This suggests that nanoislands have an amorphous or polycrystalline bulk structure, at 
least when prepared following our preparation procedure. The structure of the NPs 
does not change significantly during their sliding. Visual inspection of figure [7| indicates 
that contact surfaces of both Cu and Au NPs are not perfectly ordered but exhibit 
commensurate domains separated by domain walls. 

The results obtained above can be summarized as follows: 

• a sawtooth substrate force is observed mainly for metals with the lattice constant 
a close to the lattice constant of graphene. In our simulations it is Cu and Ni |7|; 
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Figure 6. Radial distribution function obtained at different time moments for the 
Au (a) and Cu (b) nanoparticles containing 29000 atoms. Plots for the bulk state are 
obtained using the same EAM potential. 



a 





Figure 7. Side (a, c) and bottom (b, d) views of Cu (a, b) and Au (c, d) nanoparticles 
containing 25000 atoms. Clusters of surface atoms with local order are outlined by red 
contours. 



• at low sliding velocities the sawtooth substrate force appears also for metals with 
larger discrepancy between their value of a and the graphene lattice constant. In 
our case this is true for Au; 

• RDFs indicate the absence of the long-range order in the bulk for all the metals 
prepared using our cooling procedure; 
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• visual inspection of the surfaces (see figure [?]) confirms the fact that there is no long- 
range order of the surface atoms of the NP. However, ordered domains of atoms 
can be identified on the surface of NPs. 

To explain the results of our computer simulations, we propose a "patch" model. It 
is based on the assumption of local commensur ability of some regions of the interface. As 
was mentioned above, there are domains with a local atomic order on the surface of a NP 
(see figure [7|. The domains can form commensurable regions ("ordered patches") at the 
metal-carbon interface due to some matching of metal atoms with the graphene lattice. 
The commensurate regions will effectively increase the amplitude of the corrugated 
substrate potential experienced by the nanoparticles, which will increase the friction. 
Such domains are formed more easily when the lattice constants of the NPs and the 
substrate are similar, and also for elastically |18| or plastically softer materials. We 
suggest that atomic stick-slip processes involving the "ordered patches" make the major 
contribution to the net friction force of nanoislands [7,19,20 . In general, the patches are 



randomly distributed over the interface and separated by domain walls, and are quite 
dynamic during the motion of a NP. Every patch can be considered as a locked atom and 
may be the reason for the atomic oscillations in metal. To observe a sawtooth shape 
of Fg and a higher net friction force Ff = (Fg)? all the patches must act coherently. 
This situation is more probable for Cu due to the proximity of its value of a to the 
graphene lattice constant. Such an amplification can also occur for Au and other metals 
with larger a, as was observed in the current work, when the velocity of the NP is low. 
In this case the surface atoms have time to relax into the substrate potential energy 
landscape and form pinned "patches". The extent to which such pinned domains can 
form depends not only on the lattice constants but also on the elastic or plastic properties 
of the solids and on the strength of the corrugated substrate potential ||2Tl[22l[23ll24l . 



4. Conclusions 



We have performed classical MD simulations of friction of Cu and Au NPs adsorbed 
on graphene, and sheared with a constant external force. The results obtained in this 
work confirm the earlier conclusion |7| about a significant role of the short-range order 
of atoms located on the surface of the NP. Ordered domains of atoms can be the source 
of the local match of the rubbing surfaces, leading to a sawtooth type dependence of 
the substrate force on the coordinate of the NP's center of mass. It is shown that for 
Cu, especially for small NPs, sawtooth substrate force is clearly manifested. For Au, 
which has larger a compared to Cu, sawtooth substrate force can be observed only at 
small sliding speed. The Fs{t) exhibits "fine structure", i. e. larger peaks are split into 
a set of high-frequency peaks. 

As in the previous work [t], the dependence of friction force on contact area of 
Au NPs is linear in all studied cases. The functional form of these dependencies is 
not infiuenced by changing the value of the applied force, indicating a weak velocity 
dependency, which is consistent with the experimental results (Tl[2ll3l[5] . 
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A ''patch" model is proposed in order to explain the observed behavior. Application 
of the model to different systems will be presented elsewhere. 
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